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Abstract: This paper continues the investigation of the modified Sakai-Sugimoto 
model proposed in arXiv:0708.3233. Here we discuss in detail numerical solutions 

o. 

ultraviolet cut-off turns out to be essential because the numerical solutions tend to 
rapidly diverge from the desired asymptotic solutions, beyond a sufficiently large 
value of the holographic coordinate. The required cut-off is determined by the 
non-normalizable part of the tachyon and is parametrically far smaller than that 
dictated by consistency of a description in terms of 10-dimensional bulk gravity. 



to the classical equations for the brane profile and the tachyon condensate. An 



In arXiv:0708.3233 we had argued that the solution in which the tachyon field 
goes to infinity at the point where the brane and antibrane meet has only one free 
parameter, which may be taken to be the asymptotic brane-antibrane separation. 
Here we present numerical evidence in favour of this observation. We also present 
evidence that the non-normalizable part of the asymptotic tachyon solution, which 
is identified with quark mass in the QCD-like boundary theory, is determined by this 
parameter. We show that the normalizable part of the asymptotic tachyon solution 
determines the quark condensate, but this requires holographic renormalization 
of the on-shell boundary brane action because of the presence of infinite cut-off 
dependent terms. Our renormalization scheme gives an exponential dependence on 
the cut-off to the quark mass. We also discuss meson spectra in detail and show 
that the pion mass is nonzero and satisfies the Gell-Mann— Oakes— Renner relation 
when a small quark mass is switched on. 
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1. Introduction 

The model of Sakai and Sugimoto (SS) |J has been very successful in reproducing 
many of the qualitative features of non-abelian chiral symmetry breaking in QCD. In 
this model, the 'colour' Yang-Mills fields are provided by the massless open string 
fluctuations of a stack of a large number N c of D4-branes, which are extended 
along the four space-time directions and in addition wrap a circle [[|. In the 
strong coupling limit, this stack of D4-branes has a dual description in terms of 
a classical gravity theory. Flavour degrees of freedom are introduced in the probe 
approximation as fermionic open string fluctuations between the colour branes and 
an additional set of 'flavour' branes @, f|, |5|, || [?J f§, which are provided by pairs of 
D8 and D8-branes. In this setting, chiral symmetry breaking has a nice geometrical 
picture. In the ultraviolet, chiral symmetry arises on flavour D8-branes and D8- 
branes, which are located at well-separated points on the circle, while they are 
extended along the remaining eight spatial directions, including the holographic 
radial direction. Chiral symmetry breaking in the infrared is signaled by a smooth 
joining of the flavour branes and antibranes at some point in the bulk. 

Despite its many qualitative and some quantitative successes Jj], || |H], [TT], [12], [13 



|14| , |j~5[, |i"6| , [lTj, the SS model has some deficiencies: (i) It does not have parameters 
associated with quark mass and the chiral condensate. On the other hand, there is 
a parameter, the asymptotic separation between the flavour branes and antibranes, 
which, within the SS scenario, finds no counterpart in QCD; (ii) The SS model also 
ignores the open string tachyon between flavour D8-brane and D8-brane, which may 
be reasonable in the ultraviolet where the branes and antibranes are well separated, 
but is not so at the place in the bulk where the branes join. It is often argued that 
in the curved background of the wrapped D4-branes, the geometry forces flavour 
branes to join in the interior. While this is true of flavour branes and antibranes 
that are well-separated asymptotically (separation of the order of the antipodal 
distance), it cannot be the reason when the separation is small and the branes and 
antibranes meet far away from the central region. For small separation, the effective 
radius of the direction on which the D4-branes are wrapped is very large and so 
one would expect tachyon condensation to be the primary reason for branes and 
antibranes meeting, as in the extremal D4-brane metric. Since the tachyon field 
takes an infinitely large value in the true ground state 1 , the perturbative stability 
argument given in JIJ , valid for small fluctuations of the tachyon field near the local 
minimum at the origin, does not apply. 

It has recently been suggested in |L9|, |2(], |21| that tachyon condensation on 
a brane-antibrane system describes the physics of chiral symmetry breaking in 
a better and more complete way. If the brane and antibrane are well-separated 



21 1 then one also retains the nice geometric picture of the SS model for non- 



1 For a recent review of this subject, see |TJ 



abelian chiral symmetry breaking. The purpose of the present work is to complete 



the investigations started in |y]. Here we give detailed numerical solutions to 
the classical equations for the brane profile and the tachyon. We show that the 
solution in which the tachyon diverges at the point in the bulk where the brane 
and antibrane meet has only one free parameter, which may be taken to be the 
asymptotic separation between the flavour brane and the antibrane. We present 
numerical evidence that the non-normalizable part of the asymptotic tachyon 
solution is determined by this parameter. Thus, by the usual dictionary of AdS/CFT 



[[24 |23| , |24| , pql , this parameter determines quark mass in the boundary theory 
[26| , |16 |. The parameter for the asymptotic brane- antibrane separation is present 
in the SS model also, but in that setting it cannot be explained as a parameter in 
QCD. Thus this parameter, which seems mysterious in the SS setting, finds a natural 
explanation in our model. The presence of a non-normalizable part in the tachyon 
solution necessitates introduction of an ultraviolet cut-off. This is because in this 
case the numerical solutions tend to rapidly diverge from the desired asymptotic 
solutions, beyond a sufficiently large value of the radial coordinate, determined 
by the magnitude of the non-normalizable part. This cut-off is parametrically 
far smaller than the cut-off of order iV 4//3 expected because of the breakdown of 
description in terms of a 10-dimensional gravity theory. Removing the cut-off, 
therefore, necessarily involves tuning the non-normalizable part to zero. We discuss 
how this should be done appropriately. We also discuss the chiral condensate and 
its determination by the normalizable part of the asymptotic tachyon solution. 
This determination is subtle for two reasons. One is the fact that the space-time 
independent classical solutions are described by a single parameter and hence the 
non-normalizable part of the tachyon cannot be varied independent of the other 
parameters. The resolution of this issue requires us to consider more general 
solutions by incorporating space-time dependence. But for this one has to go 
beyond the expansion in small space-time dependent fluctuations around space-time 
independent solutions, basically because this expansion is singular for the tachyon 
solution in the infrared. An exact space-time dependent action is needed, which we 
derive. The other subtlety has to do with the necessity of an ultraviolet cut-off. To 
extract cut-off independent physics, we add counter-terms to the D8-brane action to 
remove terms in the boundary action which are divergent as the cut-off is formally 
allowed to go to infinity. With an appropriate choice of the counter-terms we get 
a finite value for the chiral condensate. Finally, we discuss meson spectra in detail 
and show that the pion mass is nonzero in the presence of a non-normalizable part 
of the tachyon and that it satisfies the Gell-Mann— Oakes— Renner (GOR) relation 
when quark mass is small. 

The organization of this paper is as follows. In the next section we will briefly 
review the essential features of the modified SS model with the tachyon present. This 
section also includes a more detailed discussion of the cut-off and its implications 



than given in [21]. In section || we describe in detail the numerical solutions for 
the brane profile and the tachyon. This section also contains a discussion of the 
parameters of the solutions and their determination in terms of a single parameter, 
namely, the asymptotic brane-antibrane separation. In section [| we discuss the 
subtleties involved in deriving an expression for the chiral condensate in terms of 
the parameters of the solutions. We derive the exact 5-dimensional action in which 
the tachyon and brane-antibrane separation fields have dependence on space-time 
as well as the holographic coordinate and discuss solutions to the equations derived 
from this action. We also discuss the counter-terms required to make the chiral 
condensate finite as the cut-off is formally removed to infinity. In section |5| we 
analyse small fluctuations around the classical solution for the meson spectra. We 
show that the existence of a massless pion is guaranteed if the non-normalizable 
part of the tachyon solution vanishes. For a non vanishing non-normalizable part 
of the tachyon solution, we obtain an expression for the pion mass and derive the 
GOR relation for it. We end with a summary and discussion in section ^|. The 
Appendices contain details of some calculations. 

As this work was nearing completion, the works |28j and ^9[ appeared which 



have discussed the problem of quark mass in SS model using different methods. 



2. Modified Sakai-Sugimoto model with tachyon 

The Yang-Mills part of the SS model is provided by the near horizon limit of a set 
of N c overlapping ZM-branes, filling the (3 + 1) -dimensional space-time directions 
x^ (/i = 1, 2, 3 and 0) and wrapping a circle in the x 4 direction of radius Rk- Anti 
periodic boundary condition for fermions on this circle gives masses to all fermions 
at the tree level (and scalars at one-loop level) and breaks all supersymmetries. At 
low energies compared to Ij 1 , the theory on the D4-branes is (4 + l)-dimensional 
pure Yang-Mills with 't Hooft coupling A5 = (2ir) 2 g s l s N c , of length dimension. At 
energies lower than the Kaluza-Klein mass scale, R^ 1 , this reduces to pure Yang- 
Mills in (3 + 1) dimensions. This is true in the weak coupling regime, A5 << Rk, 
in which the dimensionally transmuted scale developed in the effective Yang-Mills 
theory in (3 + 1) dimensions is much smaller than the Kaluza-Klein mass scale, 
which is the high energy cut-off for the effective theory. In the strong coupling 
regime, A5 >> Rk, in which the dual gravity description is reliable, these two scales 
are similar. Therefore in this regime there is no separation between the masses of 
glueballs and Kaluza-Klein states. This is one of the reasons why the gravity regime 
does not describe real QCD, but the belief is that qualitative features of QCD like 
confinement and chiral symmetry breaking, which are easy to study in the strong 
coupling regime using dual geometry, survive tuning of the dimensionless parameter 
A5/ Rk to low values. 



Flavours are introduced in this setting by placing a stack of Nf overlapping 
D8-branes at the point x\ and Nf D8-branes at the point x\ on the thermal circle. 
Massless open strings between D4-branes and D8-branes, which are confined to the 
(3 + l)-dimensional space-time intersection of the branes, provide Nf left-handed 
flavours. Similarly, massless open strings between D4-branes and D8-branes provide 
an equal number of right-handed flavours, leading to a local U(Nf) L x U(Nf) R chiral 
gauge symmetry on the flavour D8 and Z)8-branes. This chiral gauge symmetry is 
seen in the boundary theory as a global chiral symmetry. 

In the large N c and strong coupling limit, the appropriate description of the 
wrapped D4-branes is given by the dual background geometry. This background 
solution can be obtained from the Euclidean type IIA sugra solution for non- 
extremal D4-branes by a wick rotation of one of the four noncompact directions 
which the D4-branes fill, in addition to wrapping the compact (temperature) 
direction. In the near horizon limit, it is given by B [27 



ds 2 = (^) {^dx'dx" + /({/) (<fa 4 ) 2 ) + (JpJ (j^- + U 2 d(l: 

where rj^u = diag(— 1, +1, +1, +1) and Uk is a constant parameter of the solution 2 . 
R is related to the 5-d Yang-Mills coupling, A 5 , which is kept fixed in the decoupling 
limit, by R 3 = ^-. Also, dVt 4 , e 4 and V4 = 87r 2 /3 are respectively the line element, 
the volume form and the volume of a unit S 4 . 

The above metric has a conical singularity at U = Uk in the U — x A subspace 
which can be avoided only if x A has a specific periodicity. This condition relates the 
radius of the circle in the x A direction to the parameters of the background by 

*-§(£)* 

For A5 >> Rk the curvature is small everywhere and so the approximation to a 



classical gravity background is reliable. As discussed in |27J, at very large values of 
U, the string coupling becomes large and one has to lift the background over to the 
11-dimensional M-theory description. 

2.1 Brane-antibrane pair with tachyon 

The effective field theory describing the dynamics of a brane-antibrane pair in a 
background geometry 3 with the tachyon included has been discussed in pTl [3i2]j. 



2 Note that U has dimension of length and is related to the energy scale U, which is kept fixed 
in the decoupling limit, by U = Ua' . 

3 For simplicity, we will discuss the case of a single flavour, namely one brane-antibrane pair. 
Generalization to the multi-flavour case can be done using the symmmctrizcd trace prescription 
of (30 . 



The simplest case occurs when the brane and antibrane are on top of each other since 
in this case all the transverse scalars are set to zero. This is the situation considered 
in |T!| . However, in this configuration one loses the nice geometrical picture of chiral 



symmetry breaking of the SS model. The geometrical picture is retained in the case 
considered in [f2(], [H[] where the brane and antibrane are separated in the compact 
x 4 direction. This requires construction of an effective tachyon action on a brane- 
antibrane pair, taking into account the transverse scalars. Such an effective action 
with the brane and antibrane separated along a noncompact direction has been 



proposed in [31, p2| . A generalization of this action to the case when the brane 



and antibrane are separated along a periodic direction is not known. However, for 



small separation l(U) compared to the radius Rk of the circle, the action in [[32 
should provide a reasonable approximation to the compact case. In the following 
we will assume this to be so. Then, the effective low-energy tachyon action for a 
D8 and _D8-brane pair for 1{U) « Rk is given, in the above background, by 5 

S = - f d 9 a V(T, Z)e~* (s/^det A L + ^-det A R ) , 
(Ai) ab = (g MN - lJ—g Mmi A d a xfd b xf + 27ra'F l ab + ^-(2na'(D a T(D b T)* + {D a r)*D b r) 

+il(g ai + d a xfg u )(r(D b T)* - r*D b r) + i/(r(D a r)* - T*D a r)(g ib - d b xfg u )J, 

(2.3) 

where 

Q = 1 + ^7^4, D a r = d a r - %{A La - A Rta )r, V(T, I) = g s V(T)^Q. (2.4) 

T = \r\, i = L,R and we have used the fact that the background does not depend on 
x 4 . The complete action also includes terms involving Chern-Simons (CS) couplings 
of the gauge fields and the tachyon to the RR background sourced by the ZM-branes. 
These will not be needed in the following analysis and hence have not been included 
here. 

The potential V(r) depends only on the modulus T of the complex tachyon r. 
It is believed that V(r) satisfies the following general properties JT8| ]: 

• V(T) has a maximum at T = and a minimum at T = oo where it vanishes. 



4 Also see [||. 

5 Strictly speaking, this action is valid only when the brane and antibrane are separated along a 
noncompact direction. However, as we shall see later, a posteriori justification for using this action 
is provided by the classical solutions for the brane-antibrane profile. In these solutions, for small 
asymptotic separation, the brane and antibrane meet far away from the central region. In this 
case, to a good approximation, the factor f(U) in the background metric can be set to identity, 
which is equivalent to setting the radius Rk to infinity. 



The normalization of V(T) is fixed by the requirement that the vortex 
solution on the brane-antibrane system should produce the correct relation 
between Dp and D(p — 2)-brane tensions. In the present case this means 
V(0) = T 8 = l/(2vr) 8 l 9 s g s , the £>8-brane tension. 

In flat space, the expansion of V(T) around T = up to terms quadratic in 
T gives rise to a tachyon with mass-squared equal to — l/2a'. 



There are several proposals for V(T) which satisfy these requirements [ 0| , 
although no rigorous derivation exists. Examples are (i) the potential used in 
P^ , [53 , pm for calculation of decay of unstable D-branes in two-dimensional string 
theory 

V(T) = T 8 sechv^T; (2.5) 

and (ii) the potential obtained using boundary string field theory computation 
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V(T) = % e"i T . (2.6) 

Both these potentials satisfy the properties listed above. Note that the asymptotic 
form of the potential in (|2.5|) for large T is ~ e~^ T . The linear growth of the 
exponent with T should be contrasted with the quadratic growth for the potential 
in ( |2.6|) . This difference will turn out to be important for the background tachyon 
solutions, which are discussed next. 

We end this subsection with the following observation. It can be easily seen 
that in the decoupling limit all factors of a' scale out of the entire action, without 
requiring any scaling of the transverse scalar I or the tachyon r. In fact, the entire 
action can be rewritten in terms of A5 and U, quantities that are kept fixed in the 
scaling limit. Henceforth, we will use the convention 2na' = 1. 

2.2 Classical equations for brane profile and tachyon 

We will now look for an appropriate classical solution of the brane-antibrane-tachyon 
system. Let us set the gauge fields and all but the derivatives with respect to U 
of T and x\ to zero. Moreover, we choose x\ = 1/2 and x A R = —1/2 so that the 
separation between the brane and antibrane is I. In this case, in the static gauge 
the action (|2.3| ) simplifies to 6 

S = -V,j 'd 4 x j 'dUV(T)(jj) f/ 4 (v^D^+v^;), (2.7) 

where Dl,t = Dr^t = Dt and 

Dt = f{U)- 1 0Q 3 2 + /([/) 0Q 3 2 ^ + T\Uf + T{U)\Uf. (2.8) 



D The CS term in the action does not contribute for such configurations. 



It is convenient to remove the dependence on R (except for an overall factor in the 
action) through a redefinition of variables, 



U = u/R 3 , l(U) = R 3 h(u), U k = u k /R 3 . 
In terms of the new variables, we get 

S = -Vyr 9 f d 4 x f du u l3/A V(T) Lfd^+^fd^ 

where 



d 



L,T — dnT 



d T = f(u)- 1 u^ + f(u)u^ )2 - - 



+ T\uY + T{uYh{u) 



(2.9) 



(2.10) 
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with f(u) = (1 — u 3 k /u 3 ). 

The effective potential for the tachyon can be obtained from this action by 
setting T' = h' = 0. It is 



V eS (T,l) ~ sechv/^zV 1 + w 3 /2T 2 /i 2 



(2.12) 



In Figure |^ we have plotted V e s as a function of T for various values of u. We see 
that a perurbatively stable minimum at T = for large values of u turns into an 
unstable maximum at a sufficiently small value of u. This is true for any fixed, non- 
zero value of h. Moreover, the value of u at which there is an unstable maximum 
at T = increases as h decreases. 




Figure 1: The effective potential V c s as a function of T for different values of u for a 
fixed non-zero value of h. 

The equations of motion obtained from the action ( |2.10| ) are 



13 
U 4 

\fd~i 



T'(u) 



13 r 
U 4 



'dn 



u 4 f(u) 
\fdh 4 



u 2 h'(u) 
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U 4 



'dn 



T(u)h(u) 2 + 

T(ufh(u) - 



V'(T) 



V(T) 

V'(T) f(u 
V{T) 4 



(d T -T'(u) 



(2.13) 



■u^h'MT'iu) 



I2.H) 



Note that the 'prime' on V(T) denotes a derivative w.r.t. its argument T and not 
a derivative w.r.t. u. 

This is a complicated set of coupled nonlinear differential equations which can 
be solved completely only numerically. To get some insight into the kind of solutions 



that are possible, however, we had analysed these equations in [|1| for large u and 
for u near the brane-antibrane joining point in the bulk. For these values of u the 
equations simplify and can be treated analytically. For the sake of completeness, 
we will summarize the results of this analysis here before proceeding to describe 
numerical solutions to these equations. As in the case without the tachyon, we 
are looking for solutions in which the brane and antibrane have a given asymptotic 
separation ho, i.e. h(u) — > ho as u — > oo, and they join at some interior point in the 
bulk, i.e. h(u) — ► at u = u > Uk 7 . Moreover, we want the tachyon (i) to vanish 
as u — > oo so that the chiral symmetry is intact in the ultraviolet region and (ii) to 
go to infinity as u approaches uq so that the QCD chiral anomalies are reproduced 
correctly |fL9|. 



2.3 Solution for large u 

Here we seek a solution in which h(u) approaches a constant ho and T becomes 
small as u —* oo. For small T one can approximate V'/V ~ — ttT 8 . If T and 
h' go to zero sufficiently fast as u — > oo such that to the leading order one might 
approximate cLt ~ w~ 3 / 2 , then ( |2.13| ) can be approximated to 

(V T'(u))' = h 2 m 4 T. (2.15) 

The general solution of this equation is 



u 



T(u) = -^{T + e- hoU + T„e hoU ). (2.16) 



In writing this solution we have ignored a higher order term in 1/u for consistency 
with other terms in equation fl2.13|) that we have neglected at large u. We will 



discuss consistency of this solution below. Let us first discuss the solution for h{u). 
The fact that the tachyon takes small values for large u makes it irrelevant 
for the leading asymptotic behaviour of h, which can be extracted from ( |2.14| ) by 
setting the r.h.s. to zero. The resulting equation is 



u -ti{u)) =0, (2.17) 



which has the solution 



h(u) = h - h lU - y l\ (2.18) 



7 The inequality results from the lower bound on u. 



'This follows from the general properties of the potential discussed in section 2.1 



Here hi is restricted to positive values so that the branes come together in the bulk. 
For SS model without the tachyon, h± = |wo/o > where f = f(u Q ), u being the 
value of u where the branes meet in the bulk. 



It is easy to convince oneself that the only solution to equations ( 2.13 ) and 



( |2.14 ) in which T vanishes asymptotically and h goes to a constant is ( 2.16 ) with 



T_ = 0. In particular, for example, these equations have no solutions in which T 
vanishes asymptotically as a power law. 

2.4 Quark mass and the ultraviolet cut-off 



In the tachyon solution ( [2.16| ), the exponentially falling part satisfies the 



approximations under which ( |2.15| ) was derived for any large value of u. The 



exponentially rising part will, however, eventually become large and cannot be self- 
consistently used. This is because for sufficiently large u, there is no consistent 
solution for T which grows exponentially or even as a power-law to the original 
equations (|2.13|) and ( [2.14|) , if we impose the restriction that h{u) should go to a 



constant asymptotically. This puts a restriction on the value of u beyond which the 
generic solution ( |2.16| ) cannot be used. The most restrictive condition comes from 
the approximation d? ~ u~ 3 / 2 . This requires the maximum value, w max , to satisfy 
the condition 



u 



5/2 
max 



rpZ -ao« maI _|_ rpZ ZhoU mBX ^ ^ IllaA /O 1 Q\ 



2K 



For generic values of \T±\ and ho, this inequality determines a range of values of 



■Umax for which the solution (|2.16|) can be trusted. The value T_ = is special since 
in this case there is no upper limit on M max , except the cut-off that comes from the 
fact that the 10-dimensional description of the background geometry breaks down 



,4/3 



beyond some very large value (~ N c ) of u. However, as is clear from ( |2.19| ) 



for nonzero |T_| one needs to choose a much smaller value of u max . Numerical 
calculations reported in the next section bear out this expectation. 

It is important to emphasize that the ultraviolet cut-off we are talking about 
here does not merely play the usual role of a cut-off needed in any example of 
AdS/CFT with a non-normalizable part present in a solution to the bulk equations. 
The point is that there is no growing solution to the tachyon equation in the 
ultraviolet which is consistent with a brane profile that goes to a finite asymptotic 
brane-antibrane separation. This constraint limits the value of u up to which the 
asymptotic solutions, ( |2.16|) and ( |2.18|) , can be trusted. 



One way to think about the inequality (|2.19|) is the following. Suppose for 
given values of \T±\ we have chosen the largest value of ii max consistent with (|2.19|) . 



Increasing w max further would then be possible only if \T_ | is decreased appropriately, 
while \T + \ can be kept fixed, as M max is increased. To be concrete, let us keep \T + \ 
and|T_|e ft ° Umax fixed as w max is increased. The process of "removing the cut-off" 
can then be understood as increasing u max and simultaneous decreasing |T_| while 



keeping \T + \ and the combination \T_\e lQUmax fixed. In this process, at some point 
|T + |e~ h ° Mmax would become much smaller than |7"_|e /l ° Umax . As we shall see in the 
next section, however, limitations due to numerical accuracy prevent us from tuning 
\T_\ to very small values, or equivalently tuning w max to be very large. Thus we are 
numerically restricted to rather small values of M max - For values of u larger than 
■u max , the inequality fl2.19| ) breaks down and consequently the asymptotic solution 



( [2.16| ) is not applicable. Clear evidence for this breakdown is seen in the numerical 
calculations reported in the next section. 

It is natural to associate T_ with the quark mass since this parameter comes 
with the growing solution. Evidence for this will be given in section [5] where we 
will show that for a small nonzero value of this parameter, the pion mass is nonzero 
and proportional to it. It is also natural to associate T + with the chiral condensate 
because it comes with the normalizable solution. It turns out that this association 
too is consistent, though this part of the story is somewhat more complicated, as 
we shall see in section f|. 

It is interesting to mention here that keeping the combination |j , _| e /loUmax = p 
fixed as the cut-off becomes large implies an exponential dependence of |T_| on 
the "Umax, i-e. |T_| = pe _/l ° Mmax . A similar dependence of the quark mass on the 
cut-off has been observed in p8| [29[, though the methods used for computing quark 



mass in these works are quite different from ours. In [^] the cut-off arises from the 



location of a D6-brane, which is additionally present in that model, thereby giving 
a physical meaning to the cut-off. 

2.5 Solution for u ~ u 

Here we are looking for a solution in which h — ► and T — ► oo as u — > uq. Let us 
assume a power law ansatz, namely 

h(u) ~ [u - u ) a } T(u) ~ [u - uo)' 13 . (2.20) 

For a smooth joining of the brane and antibrane at uq, the derivative of h must 
diverge at this point, which is ensured if a < 1. Since for this ansatz T' 2 is the 
largest quantity for u — > uq, we can approximate (It ~ T'{u) 2 . We will also need 
the asymptotic form of the potential V(T) for large T, which depends on the specific 
potential being used. From the asymptotic form of the potential in Q2.5|), we get 

V'(T)/V(T) v/?, while for the potential in (p|), we get V'(T)/V(T) ttT. 

Putting all this in ( |2.13 ) and ( p.!4j ), it is easy to verify that these equations cannot 



be satisfied by the ansatz ( (4.9|) for the potential fl2.6|) . They are, however, satisfied 
for the potential in (|2.5|) . In fact, in this case the powers as well as the coefficients 
all get fixed 9 : 



h(u) = J-^j-u 3/ \u - uo) 1 / 2 + • • • , (2.21) 



3 In p0| the power of (u — ug) with which the brane-anibrane separation falls-off in the bulk has 



been left undetermined. This power is actually determined by (2.13) and (2.14), as can be easily 



nu) = y^f ul /2 (u-u )- 2 + ---, (2.22) 

An important feature of the above solution is that it depends only on a single 
parameter, namely the value of uq. We have checked that this feature persists 
in the next few higher orders in a power series expansion in (u — Uq). This is 
in sharp contrast to the asymptotic solution ( [2.16| ), (|2.18|) which depends on all 
the four expected parameters, T + , T_, h , h\. This reduction in the number of 
parameters is similar to what happens in the SS model where the solution for u ~ Uq 
depends only on one parameter, although the asymptotic solution depends on two 
parameters. In the present case the reduction in the number of parameters is even 
more severe; the solution for u ~ Uq matches with only a one-parameter subspace 
of the four-parameter space of asymptotic solutions. As we will discuss later, this 
one-parameter freedom of the classical solution turns out to be analogous to the 
freedom to add a bare quark mass in QCD. 

For completeness, we note that there exists another solution in which T does not 
diverge but goes to a nonzero constant as u — > Uq. In this case we can approximate 
dr ~ f(u)u 3 ^ 2 h'(u) /4. Substituting in ( |2.13| ) we see that the l.h.s. diverges as 
(u — Uo)~ a - The first term on the r.h.s. vanishes as a positive power, but the second 
term diverges as (u — Wo) Q_1 , since a < 1. For consistency we must have a = 1/2. 
The resulting solution 

h(u) = -(/ (5/o + ^))- 1/2 {u - no) 1 / 2 + ■ ■ • , (2.23) 



also satisfies (|2.14j) . Note that no special condition was required for the tachyon 



potential to get this solution; this solution exists for any potential. 
3. Numerical solutions 



The equations ( |2.13| ), ( |2.14| ) cannot be solved analytically. One needs to use 



numerical tools to get a solution. We have made use of mathematica for this. 
Also, for numerical calculations we have chosen the potential (|2.5| ), since there is no 
diverging solution for T{u) for u ~ u for the potential (|2.6|) , as discussed above. 

The numerical calculations are easier to do if we start from the u = Uq end 
and evolve towards the large u end. This avoids the fine-tuning one would have 
to do if one were to start from large values of u, where the general solution has 
four parameters, and end on a one-parameter subspace for u ~ Uq. We must 



checked by consistently expanding these equations on both sides and going beyond the leading 
order in powers of (u — uq)- We have also verified this power by numerical calculations reported 
in the next section. 



also satisfy the requirement of working in the parameter region of the background 
geometry corresponding to the strong coupling. In addition, we need to ensure 
that the asymptotic separation between flavour branes and antibranes is small 
compared to the radius of the x 4 circle. Mathematically, these requirements are 



A 5 = 87r 2 i? d > 2-KRk and l < nR k . Using (fZ 
Then, these requirements become 



36vr 2 



•C u k < 



and ( 

4tt 2 



!>/'n 



one gets R = ^Rky^Uk- 
Throughout our numerical 



calculations we will work with u k = 1, which satisfies the first condition easily, while 
it requires from the second that ho <C ^. This condition is also easily satisfied by 
choosing w ~^> u k = 1 10 . For such values of uo, f(u) ~ 1 for all u> u . 

The boundary conditions are imposed using ( p.21| ), ( |2.22|) at a point u — U\ 
which we choose as close to uq as allowed by numerics. Generally we were able to 
reduce {u\ — uq) down to about 0.1 percent of the value of uq. Starting from the 
values of T(ui), T'(ui), h{u\) and h'{u\) obtained from ( |2.21| ), ( P-22J ) at u — u±, 
the system was allowed to evolve to larger values of u. Figure ^| shows an example 
for mq = 12.7. Solutions for both h(u) and T(u) are shown. 
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Figure 2: The brane profile and the tachyon solution for uq = 12.7. 
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3.1 Verification of the UV and IR analytic solutions 

From the numerical solutions one can verify that h(u) and T(u) are given by 
the forms ( |2.21| ), fl2.22| ), for u ~ uq. Figure ^ shows the impressive fits between 
the numerical data and the analytical expectations for the powers of (u — uq) 
for h(u) and T(u). We have plotted h{u)/h'{u) and T(u)/T'(u), calculated from 
the numerical solutions, as functions of u. The numerical data are plotted in 
dashed lines while the theoretical solutions are plotted in solid lines. As one 
can see, these graphs are linear at the IR end and their slopes turn out to be 
close to the expected values 0.5 and —2 respectively. In fact, the numerical and 
the theoretical curves entirely overlap in the IR region of u, as shown in Figure 
H. At the other end also, namely for large u, one can verify that the numerical 



10 As we shall see below, the asymptotic separation decreases with increasing value of uo, as is 
the case for the SS model. 
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Figure 3: Numerical verification of exponents in the IR behaviour of brane profile and 
tachyon. The fits give the two exponents respectively to be 0.50 and —2.07 for u = 13.1. 



solutions have the analytic forms ( |2.18|) , ( |2.16| ). The goodness of the fits of these 
analytic forms to numerical data is shown in Figure ^ where again the two curves 
overlap in the asymptotic region of u. The fits yield values of the four parameters: 
h = 0.224, hi = -16068, T + = 29194.5, T_ = -1.25 x 10" 4 for u = 13.1. 
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Figure 4: Numerical verification of the asymptotic form of the brane profile and the 
tachyon. 



3.2 Behaviour of the non-normalizable part 

For T_ ^ 0, extending numerical calculations much beyond the values of u shown 
in Figure ||| meets with a difficulty. It turns out that for small Uq, T_ is positive. 
Since T_ is the coefficient of the rising exponential in T(u), for a sufficiently large 
value of u this term dominates and so T{u) begins to rise n . Eventually, T becomes 
so large that the conditions under which the asymptotic solutions fl2.18j ), ( 2.16| ) 
were obtained no longer apply. Figure | illustrates this; it shows the solutions for 
Mo = 12.7 for two different large values of u. In Figure 5(a), after falling very fast, 



11 We would like to thank Matt Headrick for a discussion on this point and some other aspects 
of our numerical calculations. 
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Figure 5: Solutions for two different large values of u. 

T rises and then falls again. Almost simultaneous with this is a rapid rise of h from 
one nearly constant value to a higher constant value. Evidently, this behaviour 
continues indefinitely with u, as can be seen in Figure 5(b) 12 . 

The value of T_ decreases with increasing uq. This can be easily deduced from 
the fact that the maximum value of u up to which the asymptotic solutions ( |2.16| ), 
( |2.18p apply, namely before the oscillations begin, increases with increasing w . 
Figure |] illustrates this by showing the solutions for increasing values of uo, close to 
where T_ is small. As one can see, increasing the value of uq by a very small amount, 
from M = 13 to «o = 13.0878, dramatically increases the threshold for oscillatory 
behaviour of T from u ~ 50 to u ~ 120! As Uq increases further, T_ decreases, 
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Figure 6: Numerical solutions for increasing values of no for positive T_. 



becomes zero and eventually negative. Since we want to interpret T_ as the bare 



12 In |2(J, the authors claim that this effect is due to sensitivity of the solutions to the boundary 
conditions at the infrared end at u = U\, which must necessarily be chosen slightly away from the 
actual value uq. We have not found any evidence for this sensitivity. On the other hand, it is 
clear that the approximation made in deriving the asymptotic solution, ( 2.1 6| ) , (2.18), must break 
down for sufficiently large u, for any non-zero value of T_ . We see convincing numerical evidence 
for this. Further evidence of this follows. 

13 We have found that T_ = 1.92 x 1(T 9 at u ~ 13.0877781. Fine-tuning u such that T_ is 
precisely zero is hard. This requires numerical methods which are beyond the scope of those used 
here. However, the trend is clear from Figure @ and Figure M. 



quark mass parameter, negative values for it are allowed. However, a large value 
for |T_| will eventually again make T large in magnitude for large enough u. So 
once again we expect that at some sufficiently large u, T will become so large that 
the conditions under which the asymptotic solutions ( |2.16|) , ( J2.18Q were obtained 
no longer apply. So, as before, one should find oscillations in T(u), which now start 
at smaller and smaller u as u grows. This is indeed seen to be the case, as is 
evident in Figure H. This happens because |T_| grows with u , beyond the value 




Figure 7: Numerical solutions for increasing values of uq for negative T_ 



at which it becomes zero. Figure [8] shows the change of T_ with m . We see that 




Figure 8: T!_ as a function of uq. 



T_ vanishes at uq ~ 13.0878 and |T_| grows on both sides away from this value. 
It is hard to understand what is special about this value of uq. One might have 
thought that the role of zero mass would be played by the antipodal configuration, 
which has uq = Uk, and is beyond our approximation. It is possible that this is an 
artifact of using the approximate action, ( |2.3|) , valid for a noncompact x A coordinate, 
although the value Uq ~ 13.0878 is fairly large and seems to be within the validity of 
our approximation. We also note that for negative T_ , negative T{u) can be avoided 
by imposing a suitable cut-off on u. As we have already discussed, the cut-off is 



in any case required to fulfil the condition ( 2.19 ) so that the asymptotic solutions 
O, O may apply. 



3.3 Behaviour of the asymptotic brane-antibrane separation 

Another interesting quantity is the asymptotic brane-antibrane separation, ho, as a 
function of uq. This quantity has been plotted in Figure ||. We see that h steadily 




Figure 9: ho as a function of uq. 

decreases through the special value u ~ 13.0878. Although we do not have an 
analytical formula for the dependence of ho on w f° r large values of the latter, the 
trend in Figure |9] seems to indicate that it decreases to zero as uq becomes large. 
Presumably the brane-antibrane pair overlap and disappear as uq goes to infinity. 
This is consistent with the trend of increasing bare quark mass for increasing values 
of Uo (far beyond uq ~ 13.0878) which we have seen in Figure || Therefore, unlike 
in the Sakai-Sugimoto model, the disappearance of the brane-antibrane pair for 
Uq = oo can be understood in the present setup as the infinite bare quark mass 
limit. 

It should be clear from the above discussion that the limit ho — > does not 
reduce to the case of overlapping D8-branes and D8-branes considered in fl9| . For 
this case, one must begin afresh with x\ = 0, I = in the action ( fOl) . However, 
the classical equation for T can be obtained from the equation Q2.13J ) by setting 
h = in it. As above, we find that solutions which are divergent in the IR depend 
on only one free parameter. For further details about the tachyon solutions in this 
case, we refer the interested reader to the Appendix HI. 

3.4 Comparison with the Sakai-Sugimoto solution 

Finally, we must ensure that the solution with the tachyon has lower energy 
compared to the SS model. The energy density in the modified model is given 
by 

/*^max 

E T = 2V 4 R 9 V(0) / du E T (u), 

Juq 



E T (u) = u 13 ' A ^S-Jw 3 / 2 + \?/ 2 h'{uf + V{uf + T(u) 2 h(u) 2 , (3.1) 
while for the SS model it is given by 

/''"■max 

E ss = 2V 4 R 9 V(0) / du E ss (u), 

E ss (u) = u 13/4 Ju~ 3 / 2 + ^ U y 2 h' ss (uf. (3.2) 

To get these expressions for energy density, we have set f(u) to unity, which is a 
good approximation for large uq. Also, in the SS model one must use the solution 
of the tachyon free equation, h' ss (u) = 2mqM~ 3 / 2 (m 8 — Mq) -1 / 2 . 

Close to Mo, in the IR, the exponentially vanishing tachyon potential suppresses 
contribution to E T compared to E ss . Since the UV solutions for the two models are 
almost identical 14 , one might argue that the energy for the modified model must 
be lower than that for the SS model. However, for u > Mo there is a competition 
between the exponentially vanishing tachyon potential and the power law increase of 
the square- root factor coming from \T'\ in the integrand Et(u) in fl3.1|) . This results 
in a local maximum in E^(u) at some value of u, which can be easily estimated 
analytically. The relevant quantity, 

e -i«r ( «-« o)- 2 (M _ Mo) - 3) 

has a maximum at m = m + (| ) 1//2 w . For small m , the position of the maximum 
is close to Mo, so in this case the argument about the IR behaviour of the integrand 
in ( p.ip is not very clean, except in the very deep IR. But since the position of the 

3/4 

maximum grows with increasing Mo as m , our argument should hold for large values 
of Mo, which is precisely where the action for the modified model can be trusted. 
However, the expression used for estimating the position of the local maximum 
breaks down if it is too far away from m . So, in practice we need to do a numerical 
calculation to see what the real story is. As we will see in the numerical plots given 
below, what really happens is that for relatively large values of Mo the integrand 
Et(u) increases rapidly at first, then slows down almost to a constant and finally 
settles into an asymptotic power law increase similar to that of the integrand Ess(u) 
for the SS model. Moreover, the place where the rapid increase begins shifts to larger 
values of u as mq increases, in accordance with the above expectation. 



We have numerically evaluated the integrals in (|3.1| ) and ( |3.2| ). Because the 



relation between mq and the asymptotic brane-antibrane separation is different in 



14 There is a caveat here. Strictly speaking this is true only when the coefficient of the non- 
normalizable term, T_, in the asymptotic tachyon solution ( pi.lq ) vanishes. As we have discussed, 
when T_ is nonzero, one must introduce a cut-off, u max , chosen carefully such that the asymptotic 
solution is satisfied. In particular, one must ensure T is positive in the region below M max - In the 
calculations reported here and earlier in this section, this is what we have done. 



the two models, a given value of uq corresponds to two different values of the latter 
and vice versa. We have chosen to do the comparison for the same value of the 
asymptotic brane-antibrane separation in the two models, but the conclusions are 
similar with the other choice as well. In Figure |l(] we have plotted numerical 
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Figure 10: h(u) and T(u) profiles for u$ = 17. For comparison, /tgs profile has also been 
plotted after adjusting the value of uq to 16.4 for it since this value of no produces the 
same asymptotic brane-antibrane separation. 




Figure 11: The energy density integrands E$$(v) and E^(u). The rapid rise of the latter 
in the IR is clearly seen. The divergence between the two curves in the asymptotic region, 



u 



> 



is due to a nonzero T_ 



solutions for h(u) and T(u) for Uq = 17 15 . For comparison with the SS model, 
we have also plotted hss after adjusting the value of Uq for it to produce the 
same value of the asymptotic brane-antibrane separation. The required value turns 
out to be m = 16.4. The corresponding energy density integrands, Ei(u) and 
Ess(u), have been plotted in Figure [11]. We can clearly see the rapid rise of 
E(u) in the IR, the subsequent flattening out and finally the power-law rise in 



15 Similar behaviour is seen for values of u > 14. Below w ~ 14, however, the energy difference 
becomes very small and even reverses sign. This may be connected with the breakdown of the 
approximate action in this region, similar to the observation of a zero quark mass at uq ~ 13.01. 



the asymptotic region. Using w max = 35.32 16 , numerical evaluation of the integrals 
gives (Et — Ess) = —300.3. Therefore, the solution with the tachyon taken into 
account corresponds to a lower energy state. 

4. The chiral condensate 

By the standard dictionary of AdS/CFT |25], [22|, ^3], |23]|, once we have identified T_ 
with the quark mass parameter, we should identify T + with the chiral condensate. 
However, it is not clear that the standard rules apply to the present case of a 
boundary theory which is not a CFT and has a scale. Moreover, the fact that 
there is no known lift of D8-brane to 11-dimensions forces an essential cut-off in the 
theory with flavours. In fact, for a non-zero value of T_, the real cut-off is much 
lower, as we have seen from numerical computations in the last section. Despite 
these difficulties, we will assume that the identification of sources in the boundary 
theory with boundary values of bulk fields holds in the theory with cut-off. 

There is an additional difficulty in the present case. As we have seen above, the 
desired solutions have only one independent parameter, which we take to be T_. The 
other three parameters, T+, ho and hi should then be considered to be functions of 
T_ . Thus, the chiral condensate cannot be computed naively by varying the on-shell 
flavour brane action w.r.t. T_, since this would also include contributions from the 
variation of the other three parameters with T_ . The one-parameter solutions that 
we have found constitute the most general class of space-time independent solutions 
with the specified boundary conditions 17 . Therefore, if we only want to make a 
variation of T_ only, we must go out of the present one-parameter class of solutions 
to more general solutions, which are space-time dependent, in addition to being 
dependent on u, and have enough parameters. These solutions to (u, x)-dependent 
equations should have the same singularities at u = u as the solutions in (|2.21| ) 
and ( p.22|) . Moreover, the asymptotic solutions should have the form of Q2.16|) and 



( [2.18|) with x-dependent coefficients. If solutions satisfying these conditions exist 
and have enough parameters, then we can make the required variation of T_ only 
and identify T + as the condensate in a coherent state formed from fluctuations of 
T and h (scalar mesons) around the ground state with broken chiral symmetry. 
Specializing to the ^-independent case, after varying the on-shell action, then, gives 
us the condensate in the vacuum state. What we, therefore, need to do is to analyse 
the ^-dependent case to see if the required solutions exist. This is what we will do 
next. 



16 



This is the value at which T{u) vanishes. The asymptotic form, (2.1 6|) , fits the numerically 



computed T(u) in the range 33 < u < u max to better than a percent with the parameter values 
h = 0.179, T+ = 28904, T_ = -0.0937. 

17 These boundary conditions are (i) vanishing tachyon and fixed brane-antibrane separation 
asymptotically and (ii) divergent tachyon and vanishing brane-antibrane separation at some point 
in the bulk. 



4.1 Action for (u, x)-dependent T and h 

The full (u, x)-dependent action for tachyon and brane-antibrane separation is given 

by _ 

S = -^ I d 4 x f duu 13/A V(T)y/ck y/det(l + K), (4.1) 

where K is the matrix with the elements 

f u~ 3 / 2 

K ^ = {o d>lhd » h + -Q- d " Td » T ' 

f u- 3/2 

fu 3 ' 2 1 
K u = — h'B h ^ T'B T 

K\ = 0. (4.2) 



To look for a generalization of the x-independent solutions for equations of motion 
derived from this action, the most obvious thing to do is to generalize the earlier 
solutions by making all parameters functions of x. In particular, this means making 
Mo, the place where the flavour brane and antibrane meet, a function of x. For 
u ~ Mo, expansion of this solution around a constant Uq is singular, since it involves 
arbitrary higher powers of l/(u — Uq). Therefore, we do not expect analysis of ( |4.1|) 
by expanding in small fluctuations around the x-independent solution to work for 
u close to Mo- This is confirmed by explicit fluctuation calculations in Appendix [B| 
We need to go beyond small fluctuations analysis of fl4.1|) and this requires us to 
get an exact expression for the determinant in terms of space-time derivatives of T 
and h. 

A direct calculation of det(l+A) is tedious, but the calculation can be simplified 
using a trick which has been described in Appendix 0, where a rather simple 
expression for the determinant has been obtained. The complete 5-dimensional 
action then reads 

S = -^ J d A x I du m 13/4 V(T) v/A^, (4.3) 

where A^ = <irA and we have defined 

A = 1 + p^dT) 2 + P 2 (dh) 2 + 2p 3 (dh.dT) + p 4 [(dT) 2 (dh) 2 - (dh.dT) 2 }. (4.4) 

The /3's are given by 

-3/2 T ,2 f f u 3 / 2 h' 2 fh'T' 

^-— {l -ir^ e^k 1 -^' ft = -ik' *=aa-# (4.5) 



As a check on the action (|4.3| ), we note that it reduces to the action ( |2.10|) if T and h 



are x-independent. Also, it correctly reproduces the action (|B.1|) which only retains 



terms that are quadratic in space-time derivatives of T and h. This latter action 
was derived independently by expanding det(l + K) in powers of K and retaining 
only the first nontrivial correction. 

The equations of motion that follow from the action ( [Of ) are rather complicated 
and have been derived in Appendix |C], ( |C.13| ) and ( |U.14| ). As we did in the x- 
independent case, we will solve these equations in the two limiting cases of large u 
and u ~ Uq. 



u 



u. 



: In this limit, h(u,x) goes to a fixed value ho(x), which is assumed to 
be a slowly varying function of x. We will also assume that T and all its derivatives 
are small in this limit. Then the equations (|C.13| ) and ( |C.14| ) can be approximated 



as 



(u 4 T'(u, x))' + (h (x)) 2 u 4 T(u, x) = 0, 

u^~h'(u, x) ) = 0. 



(4.6) 
(4.7) 

The space-time derivatives are comparatively suppressed by powers of 1/u and hence 



have been ignored. These equations are identical to (|2.15|) and (|2 . 1 7|) and so have 
solutions similar to ( 2.16Q and ( |2.1S ), but now with parameters that are functions 
of x: 



T{u,x) = \{T + {x)e~ ho{x)u + T_{x)e ho{x)u ), 
h(u,x) = h (x) — hi(x)u~ 9 ' 2 . 



(4- 



u 



Uq'. The analysis in this limit is somewhat more involved. We assume an 



ansatz similar to the solutions (|2.21| ) and (|2.22| ), but now with x-dependent Uq and 
coefficients: 



h(u,x) = po{x)(u - uq(x)) 1/2 + pi(x)(u - uq(x)Y 12 + 
T(u,x) = a (x)(u — u (x))~ 2 + <Ji(x)(u — u (x)) _1 + ■ 

As consequence of this ansatz, one can show that 

d^h = -h'[d,j,u — (u - w ) H ], 

Po 

d^T = -T'[d,u + ^(u-u ) + •••]. 
2cr 



(4.9) 

(4.10) 
(4.11) 



These relations are correct to the order shown. Putting all this in the equation 
of motion for T, ( |U.13|) , we see that this equation is satisfied to the leading order 
provided the following condition holds: 



13 

Au 



°oPo = M 3/2 ^( M 3/2 ^ W o) 



-ufd»u 



d,{u- \du Q ) 2 ) 



(4.12) 



In obtaining this we have set fo — 1. Similarly, from QC.14 ) one gets the condition 



on = ^(^ /2 H-< 3/2 (^o) 2 ). (4.13) 



If Mo is a constant independent of x, then from equations ( f4.12| ) and (|4.13|) one 
gets 



v^ 3/2 / 26 _ 3/4 

ct = -j-V > Po = W % • 4.14 

4 V ^'"o 



These reproduce the x-independent solutions in fl2.21| ) and (|2.22j ), remembering that 



we have set /o = 1. Let us now consider a small fluctuation around this constant 
solution. Linearizing the equations Q4.12 ) and ( |4.13 ) in fluctuations, we get 



oo"o(x) = — — u ' du (a;), opo(^) = /==(» + -^ M o)£>wo- (4.15) 

o V2d7T o 

Now, clearly we could choose the fluctuation 5uq(x) to be such that Spo(x) vanishes. 
Under such an infinitesimal change of u , a would change, but not p . It is this kind 
of greater freedom in independently varying the parameters of the solution that we 
have wanted. Presumably in higher orders the situation gets better because there 
are more terms in the ansatz (|4.9|) and for each coefficient there is some freedom 
because of the space-time dependence. It would be nice to analyse the higher order 
terms, but that is beyond the scope of this work. Here we will assume that the 
introduction of space-time dependence as above can give us the required freedom 
to do the calculation of the condensate as follows. 

Finally, let us compare the solution (|4.14|) , ( |4.15| ) with the solution obtained 
by the singular perturbation expansion in Appendix |B], ( |B.31| ). Expanding 



around constant u$ solution to the lowest nontrivial order in e = (u — «o) and 
comparing with ( p.28|) , we get the relations 

<p (x) = 2 Su (x), tpi(x) = —(Sa (x) + a 1 8u (x)), 

1 13 

& Q (x) = --5u Q (x), 0i (a;) = — (5p Q (x) - -p x Su (x)). (4.16) 

i po z 



These relations involve not only the leading order parameters (|4.14f) of the constant 



solution, but also the nonleading parameters 0\, p%, which are given by 

^i = 7— > Pi = -5— ■ (4.17) 

6Mq oUq 



Using (|4.14j) - (|4.17| ), one can show that the equations in (|B.31|) are satisfied. This 



equivalence is, however, only formal. As we have argued above, the method given 
in this section is the correct one to use since it does not involve a singular expansion 
in arbitrarily high powers of l/(u — «o). 



4.2 Condensate in terms of the tachyon solution 

To derive an expression for the condensate, we calculate the variation of the action 
in (|4.3| ) under a general variation of T and use the equation of motion ( |C.13|) to 
reduce it to a boundary term: 

^ = -^ [ ^x^^T'(u,x)6T(u,x)\ u=u _. (4.18) 



We have ignored terms with space-time derivatives because from now on we will be 
specializing to the x-independent case, except in the variation 5T, so these terms will 
drop out. Only the UV boundary contributes to the on-shell action; there is no IR 
contribution because the tachyon potential vanishes exponentially for the diverging 
tachyon in the IR. We are only interested in retaining the variation <5T_ (x), so we set 



5T + (x) to zero. Doing this and using (}4.8|) in (|4.18| ), we get the leading contribution 
for large w max , 

SS « 2hoV ^ {0) (T + - T_e 2 ^— ) fd 4 x 6T.(x). (4.19) 

On-shell brane actions have UV divergences which need to be removed by 
the holographic renormalization procedure 18 to get finite answers for physical 
quantities. One adds boundary counter terms to the brane action to remove the 
divergences, following a procedure described in | i3| . Our on-shell action ( |4.19| ) 



diverges as the cut-off is removed. This is because, as discussed in section pT4| , we 
are keeping T + and y_e h ° Umax fixed as the cut-off is removed and the last term in 
( |4.19|) diverges as e h ° Umax in this limit. The holographic renormalization procedure 
has been developed for examples with CFT boundary theories. Since, with the 
D8-branes present, there is no 11-dimensional description available to us, it is not 



clear that the procedure described in [43 is applicable to the present case. We will 



proceed on the assumption that this is the case. Therefore, to subtract the UV 
divergent term in ( |119j ), we will add the following counter term to the boundary 
action, 

Set = ^jP-Jd 4 x V^ h(u,x)T\u,x)\ u _ } (4.20) 

where 7 = -m^ is the determinant of the metric on the 8-dimensional boundary 
orthogonal to the slice at u = w max . Note that the counter terms must be even in 
powers of the tachyon because of gauge symmetry. Using the solution (|4.8| ) and 
retaining only the parameter T_(x), we find that the variation of the counter term 
action is 

SS ct = 2fe °y 0) (r + + r_e 2 ^) f d A x ST_(x). (4.21) 



18 For reviews, see jll], Q. 



Adding to ( 4.19 ), the divergent term drops out and we get the variation of the 
renormalized action 



SS TI 



W (o), 

W 



T, / <Px5T(x) 



(4.22) 



Note that the variation of the renormalized action is twice as large as it would have 



been if we had simply dropped the divergent term 19 in (|4.19|) . 



We are now ready to calculate an expression for the chiral condensate in terms 
of the parameters of the tachyon solution. The parameters T± are dimensionless. 
To construct a parameter of dimension mass from T_, we introduce a scale /i and 
define m q = /i|T_|. Then, identifying the chiral condensate x — < QlQr >, with 
SS renorm /fi6T^(x), we get 



X 



4/ioVkV(0) 



t+ 



(4.23) 



We see that the parameter T + determines the condensate. Figure |T2| shows a plot of 
Ti as a function of T_ for T_ ~ 0. T. seems to attain a maximum value at T_ = 



"0 



T+ 

60000 r 




Figure 12: T + as a function of T_. 



and drops off rapidly, at least for small values of \T_ 



19 In (4.19), it is inconsistent to drop the term proportional to T_ in the limit of large cut-off, 
holding T + and T_e hoUnxax fixed. In fact, it is the T_ term that dominates in the action ( 4.19 ) in 
this limit. Taking a different limit that allows one to simply drop this term creates difficulties in 



the calculation of the pion mass, see section |5.3| . Consistency with the chiral condensate calculation 
then demands that the term proportional to T + T- be dropped in the pion mass calculation since 
it is smaller than the T? term. 



5. The meson spectra 



In this section we will discuss the spectra for various low spin mesons which are 
described by the fluctuations of the flavour branes around the classical solution 



20 . The action for the fluctuations of the gauge fields can be computed from (2.3) 
Parametrizing the complex tachyon r in terms of its magnitude and phase, r = Te %e 
we get the following action, correct to second order in the fluctuations: 



AS, 



gauge 

a{u) 
b{u) 

c(u) 

d{u) 



d 4 x du 



a\u) 



+d(u) ((F;j 2 + (i^) 2 ) 
R- lb V A V(T)u l 



flU/ 

T 2 



)K + b(u)Al + c{u) {(Fp 2 + (i^) 2 ) + e(v)F%A" 

(5.1) 

(5.2) 

, T 2 / f 2 T 2 h 2 h' 2 \ 

R-^ViTW^h 1 - (l + '\* h u 3 ) , (5.3) 



fdn 



Adi 



R 3 



V 4 V(T)u 1/4 ^/d 7 



R-VV A V(T)u 7 /*-Q- 



e(u) = R~ 6 V 4 V(T)u 



13/4 



4v3t 
fT 2 hh' 



1\fdn 



(5.4) 
(5.5) 

(5.6) 



Here F^ v is the usual field strength for the vector gauge field V = (Ai + A 2 ) and 
F^ v is the field strength for the gauge-invariant combination of the axial vector field 
and the phase of the tachyon, A = (A\ — A 2 — 89) . However, 



fJ,U 



UfJ, 



d,y u - R 3 d u v 



A'" 



flU 



Upb 



(JnA u — R O u An. 



(5.7) 



The relative factor of R 3 simply reflects the change of variables (2.9). 

The gauge field V yk {x 1 u) gives rise to a tower of vector mesons while the fields 
An(x,u) and A u (x,u), which are gauge invariant, give rise to towers of axial and 
pseudoscalar mesons. Notice that the coefficients a(u), b(u) and e(u) vanish if the 
tachyon is set to zero. In the absence of the tachyon, the vector and axial vector 
mesons acquire masses because of a nonzero d(u), but there is always a massless 
"pion" 21 . The presence of the tachyon is thus essential to give a mass to the "pion" . 
Also note that with the tachyon present, the masses of the vector and axial vector 
mesons are in principle different. 



20 For a general review of mesons in gauge/gravity duals, see J44|. 

21 Strictly speaking, for the U(l) case under discussion, this pseudoscalar is the v( '. It is massless 



here because of the N c — > oo limit in which we are working. 



5.1 Vector mesons 

We will be using the gauge V u = 0. Expanding in modes, we have 

V,(x,u) = J2vt ] ( x ) W m(u), (5.8) 

III 

where {Wj^tt)} form a complete sets of basis functions. These satisfy 
orthonormality conditions which will be determined presently. The fields {Vp i (x)} 
form a tower of vector mesons in the physical (3 + l)-dimensional space-time. In 
terms of these fields, the vector part of the action ( |5.1| ) takes the form, 



A^ auge = ~fd 4 xJ2 [Ql n F^ m) F v ^ + L v mn V^V^] , (5.9) 

where F^ are the usual (3 + l)-dimensional L/(l)-invariant field strengths for the 
vector potentials {V/, }. Also, we have defined 

Q v mn = [du c(u)W m (u)W n (u), L v mn = R 6 [du d{u)W' m {u)W' n (u) . (5.10) 



mn r> 



i? 6 (d{u)Wi{u)W n {u)) + 2XlQl n ] +m~n, (5.12) 

V / au 



In addition, we choose the basis functions {W m {u)} to satisfy the eigenvalue 
equations 

-R 6 {d{u)W' m {u))' = 2\ v m c{u)W mi (5.11) 

Using these we see that 

in<n„)\\ ;j„)\\ „(//); 

du 

where, as in the previous section, du refers to boundaries in the u-direction. 

Note that a potential zero mode in the vector sector 22 can be gauged away using 
the residual symmetry of making ^-independent gauge transformations, which is still 
available after fixing the gauge V u = 0. This is because a zero mode in this sector 
can only have a single scalar degree of freedom. This follows from the requirement 
of finiteness of the the action, ( |5.9|) , which cannot be satisfied since the coefficient 
of the field strength term blows up for a zero mode. Hence its field strength must 
vanish, leaving behind only a longitudinal degree of freedom. 

For the nonzero modes we may, without loss of generality, choose 

Qmn = ^Omui (5.13) 



which, on using Q5.12J ), gives 



L mn — 2^m$mn- (5-14) 



Using (£T3|) and (|53D in Q, we get 

A5j auge = - J d'x J2 [\f^F v ^^ + \\lV^V^] . (5.15) 



22 



A zero mode is defined as a mode which has zero eigenvalue and goes to a constant at infinity. 



5.2 Axial vector and pseudoscalar mesons 

As we have already noted, A^ and A u are gauge invariant. Expanding in modes, we 
have 

A,(x,u) = Y J ^\x)P m {u), A u (x,u) = J2<P (m \x)Srn(u), (5.16) 

m m 

where {P m {u)} and {S m (u)} form complete sets of basis functions. These satisfy 
orthonormality conditions which will be determined presently. The fields {A™ (x)} 
and {(fi( m \x)} form towers of axial vector and pseudoscalar mesons in the physical 
(3 + l)-dimensional space-time. In terms of these fields, the axial- vector and 
pseudoscalar part of the action (5A) takes the form, 

AS g A augc = - fd 4 xJ2 [^ mn At0 (m ¥ n) + QLK (m)FA{n)fiV + Li n A^A^ 

m,n 

+K mn d^d^ + J mn A^d^] , (5.17) 

where F^y are the usual (3 + 1) -dimensional L/(l)-invariant field strengths for the 
axial vector potentials {A}T '}. Also, we have defined 

Qmn = I du c(u)P m (u)P n (u), 

Li n = J du (B?d{u)P' m {u)P' n {u) + (b(u) + l -R z e'{u))P m {u)P n {u)), 
Jmn = [du (e{u)P m {u)-2R i d{u)P' m {u)}Sn{u), 
K mn = du d(u)S m (u)S n (u), (5.18) 



and used the orthonormality condition in the pseudoscalar sector 

du a(u)S m (u)S n (u) = -Xfa^mn- (5.19) 

In addition, we choose the basis functions {P m (u)} to satisfy the eigenvalue 
equations 

-R & {d{u)P' m {u)^ + (b(u) + l -P}e'{u))P m {u) = 2\ic{u)P m {u). (5.20) 

Using these we see that 
1 



T A 

ran r> 



R 6 (d{u)P' m {u)P n {u)) + 2\iQi n ] +m~n } (5.21) 

V / au 

where, as before, du refers to boundaries in the u-direction. 



We note that because of the last term in ( 5.17 ), the longitudinal component of 
A™ and ( - m - ) mix. So we need to define new field variables in terms of which the 
action ( |5.17| ) is diagonal. Before we do that, let us first note that the axial vector 
potential A^(x,u) has a possible zero mode provided the corresponding (3 + 1)- 
dimensional field strength vanishes, for reasons explained in the previous subsection. 
Hence the zero mode, which we shall denote by A^ , can only have a longitudinal 
component. The zero mode is gauge- invariant and, because of its mixing with the 
pseudoscalars, plays a special role. Let us see this in some detail. 

The zero mode Ajj, is conjugate to the eigenfunction Pq(u) which satisfies the 
equation 

-R 6 (d{u)P^{u)^ + (b{u) + ^R 3 e'(u)^P (u) = 0. (5.22) 

If there is no solution to this equation, then the zero mode does not exist and we 
should proceed directly to diagonalize the action ( |5.17|) . If, however, a solution 
Po(u) to this equation exists and is such that it goes to a constant at infinity, then 
the zero mode A^ exists. Since it is purely longitudinal, for a reason identical to 
that discussed in the vector case, we make this explicit by writing it in terms of a 
pseudoscalar field, Ajj. = d^4>^\x). The terms in the action ( |5.17|) which contain 
4>^\x) can be separated out. These terms are: 

in 

The sums over the indices m, n no longer include the zero mode. Also, we have used 
L^ = L^ m = for m ^ 0, which follows from (|5.21|) using the fact that A^ = 
and the boundary terms vanish because P m (u) vanishes sufficiently fast at infinity. 
Without loss of generality, we may choose Lq = 1/2 (to get the normalization of 
the kinetic term of (p^ right). Then, we can rewrite the above as 

\dpitd»it ~IJ2 Jo m Jond^ {m) d^ n \ (5.23) 



2 p 2 



rnji 



where 7T = (0(°) + E m ^o m (m) ). 

With the zero modes explicitly separated out in this way, for the nonzero modes 
we may, without loss of generality, choose 

1 



which, on using (|5.21| ), gives 



Qin = 4<W, (5.24) 

L m n = 2^m^mn- (5.25) 

Putting fT23l) , QQ4D and (|Q5D in the action (|5T7| ), we get 

m 
+ l -8^d^ + Y,{Kmnd^ {m) d^^ + JmnA^d^)] , (5.26) 



rn.n 



where K mn = (K mn — \Jo m Jon)- The above action describes a massless particle, 
7r, besides other massive particles. The existence of this massless particle depends 
on the existence of a solution to the equation fl5.22| ), satisfying the normalization 
condition 

R 6 (d(u)P (u)P^u)) du =-. (5.27) 

Later we will see that the existence of the desired solution Pq(u) depends on the 
absence of a non-normalizable part in T(u). 

To diagonalize the action ( |5.26| ) for the massive modes, we define the new 
variables 

A M = A (rn) _ J2(Xir l J mn d^ n \ (5.28) 

n 

Putting in ( |5.26| ), we get 

A^luge = - / d 4 x [J2 (~A£</> (m V (m) + \f^F a ^» + l -X A m A^A^ 



Id^d^ + ]T K' mn d^d^ n) 



(5.29) 



where K' mn = [K mn — | Yl p (^p) 1 Jpm,Jpn)- The modes have now been decoupled. To 
get the standard action for massive pseudoscalars we may, without loss of generality, 

S6t /I 1 i 



n ~ 1 1 II L I I I II. n ~ \JI Itr \Jlti f-i 

P 



This condition can be rewritten in a more conventional form as follows. We define 
^ m (u) = Y,{K)~ 1 Pn{u)J nm + P (u)J 0m , (5.31) 

n 

and using ( 5.2 0| ) note that it satisfies the equation 

-R 6 ^d(u)iP' m (u}^+ (b(u) + ^ 3 e'(«))Vm(«) = ^e(u)S m (u) + R 3 (d(u)S m (u))' ■ 

(5.32) 
Using ( |5.31|) in ( |5.30|) , we get 

5 mn = / du(d(u)S m (u)(S n (u) + R 3 ip' n (u)) - -e(u)S m (u)tp n (u)) + m <-> n. (5.33) 

In terms of new variables defined by 

S m (u) = R 3 r]' m (u), 9 m (u) = ip m (u) + r) m (u), (5.34) 

( |5.33| ) can be written as 

du rfju) (R 6 d{u)9' n {u) - ^R 3 e(u)(6 n (u) - »j n ( u ))) + m *-> n = 5 mn . (5.35) 



Moreover, in terms of these variables the differential equation (|5.32|) reads 

-R G (d{u)9' m {u)) + (b{u) + l -R 3 e\u)) (e m {u) - r) m (u)) - ±R 3 e(u) V ' m (u) = 0, 

(5.36) 
From these two equations one can obtain the orthonormality condition 



du(R b d(u)6' m (u)6' n (u) + (b(u) + -R 6 e'(u))(9 m (u) - r] m {u)){9 n {u) - Vn (u)) 



--R 3 e(u)r]' m (u)(9 n (u) - rj n (u)) - -R 3 e(u)r]' n {u)(9 m (u) - rj m (u))) = -S mn . 



(5.37) 



Also, rewriting ( |5.19| ) in terms of the new variables, we have 



R 6 I du a(u)r]' m (u)r]' n (u) = -A^<W 



(5.38) 



Finally, ( |5.35|) and Q5.38J ) give 



R 6 a(u) V ' n (u) = \i(R 6 d(u)6' n (u) - -R 3 e(u)(6 n (u) - Vn (u)) 



(5.39) 



Equations ( |5.36| ) and ( |5.39| ) are the final form of the eigenvalue equations and ( |5.37| ) 
and ( p.38[ ) are the orthonormality conditions in the pseudoscalar sector. 

It is interesting to note from ( |5.36|) that if r\ is constant, then the variable {9—rj) 
satisfies a differential equation that is identical to the equation (|5.22|) satisfied by 
the zero mode Pq. Also, using (|5.36|) and (|5.37|) one can show that for constant i], 
(9 — i]) satisfies the normalization condition ( |5.27| ). From ( |5.39| ) it follows that if 77 is 
constant, the eigenvalue A* vanishes. Thus, the presence of a massless pseudoscalar 
can be naturally considered to be identical to the question of the existence of a 
solution to the equations ( |5.36| )- fl5T39| ) with zero eigenvalue, and so it becomes a 
part of the spectrum in the pseudoscalar tower of states. Hence, the action in this 
sector can be written in the form 



AS A 

gauge 



f d 4 x X^\\f^F a ^ v + -X^A^A^ 






(5.40) 



Note that we have dropped the field tt(x), but extended the sum over m to cover a 
possible zero mode as well. If there is a solution to the equations (|5.36|) - (|5.39|) with 
constant 7/0 and Aq = 0, then a massless pion field will reappear as the zero mode 
(jy® in the pseudoscalar tower. Otherwise, the lowest mode in this sector will be 
massive, whose mass can be computed as in the following subsection. 



5.3 Relation between pion mass and non-normalizable part of tachyon 

In this subsection we will derive a relation between the pion mass and the non- 
normalizable part of tachyon parametrized by T_ . This will give us further evidence 
for identifying the parameters T + and T_ with the chiral condensate and quark 
mass respectively. We first note that for T{u) = 0, a(u) vanishes and hence A^ 
also vanishes by (|5.39|) . However, as we will see from the following calculations, 



T(u) = is a sufficient condition, but not necessary to guarantee the presence of a 
massless pion. The necessary condition is that the non-normalizable piece in T(u) 
should be absent, i.e. T_ = 0. 

Let us assume that T(u) ^ so that a{u) ^ 0. Then, ( |5.39|) can be used to 



solve for r] m (u) in terms of ip m (u), which is related to 6 m {u) and r) m (u) by (|5.34j) . 
We get, 

V'm(u) = ^^ A d{uWJu) ~ e ^m{u)) (5.41) 

a(u) - \md{u) v 2i? d / 

Let us now denote by Aq the lowest mass eigenvalue. The corresponding 
eigenfunctions are ipo(u) and t]q(u). Assuming Aq <C a(u)/d(u) 23 , we can 
approximate the above equation for rj (u): 

V' (u) « ^y (d(u)i&(u) - ^Uv)) (5-42) 

If we know i/)o(u), then using the above in ( |5.38| ) we can compute the mass. Now, 
ipo{u) satisfies the following differential equation, which can be obtained from ( 5.36Q 
using ( |5.42| ) and the approximation Aq -C a(u)/d(u): 

-R 6 (d(u)ij)' Q (u)) + (b(u) + -R 3 e'{u))4j {u) « 0. (5.43) 



Also, using (|5.43 ) and the approximation under which it was obtained, the 



normalization condition on ipo(u) given by ( [5.371 ) can be approximated as 

R*d{uW (u)^{u)\ u=Ul ^^ 1 -. (5.44) 

These equations cannot be solved analytically in general. However, analytic 
solutions can be obtained in the IR and UV regimes. In the UV regime, for u < w max , 
we use ( |2.16| ) and ( |2.18| ) to approximate the coefficients in ( |5.43| ); we get 



,, , V A V(0) m2r N ,, N V 4 V(0) 5/2 , , 9^(0), . -3/2^2/ n 



(5.45) 



23 This approximation can be justified a posteriori by the solution because the eigenvalue Ag 
turns out to be parametrically much smaller by a factor of 1/i? 3 , see (5.56), compared to the ratio 
a(u)/d(u). 



In writing these, we have used f(u) ~ 1, which is a good approximation for large 
u. We see that we can clearly neglect e{u) compared to b{u) in ( |5.43|) , while b{u) is 



itself negligible compared to d(u). Using these approximations in ( |5.43| ) and ( |5.44j ) 
then gives 

-(«*$(«))' «0, ^^«^o(«)^)(«)l«=w«i, (5-46) 

which are solved by 

*<"> " * " 3^^) , '" 3/2 ' (5 ' 47) 

Here cq is a parameter which is related to the pion decay constant. This can be 
argued by analysing the 4-d axial current correlator and using AdS/CFT along the 
lines of [[45], |46|. Using the AdS/CFT dictionary, one can compute the axial current 
correlator from the action ( |5.17|) , evaluated on-shell, by differentiating twice with 
respect to the transverse part of the axial vector field on the UV boundary. This 
is the source which couples to the axial current on the boundary. The source arises 
from the same zero mode solution, Po{u), which we discussed in connection with 
a possible zero mode (the pion) in the longitudinal component of the axial gauge 
field. Pq(u) satisfies the equation ( p.22| ), which is identical to that satisfied by 



ipo(u), (|5.43D. However, the boundary condition now is different; it is the boundary 



condition for a source, -Po(w m ax) = 1. In addition, one imposes the condition 

R e d(u)P^u)P (u)\ u=Ul ^ « & (5.48) 

which is required to reproduce the correct zero momentum axial current correlator 
|^, ^6[. This follows from the action ( p. 17] ). Now, Pq(u) satisfies ( 5.22 ) and the 



condition ( p. 48] ) if we set Pq(u) = f^oiu). Then, requiring -Po(^max) = 1 gives 

Co = I//*. 

In the IR regime, u > uq, we use ( 2.21| ) and (|2.22| ) to approximate the 



coefficients in ( |5.43|) ; we get 



IX 



3 / 2 V±u 7/i V(T) lW^\ nrr , , v lW A uT V(T) 



b(u)fa ^ — k -^tj, d(u)m ^-rV(T), e(u).~ 

K ! 26R 3 {u-u Q ) Al v ; 32^R 9 K h K J 16^/ttR 6 (u - u ) 

(5.49) 
In writing these, we have used /(wo) ~ 1, which is a good approximation for large uq. 
Using dV{T)/du = T'(u)V'(T), we see that b(u) and R 3 e'{u) both go as (u - u )~ 4 
in this regime. However, the coefficient of the latter is suppressed by a relative 
factor of Uq , so for large u we may neglect it compared to b(u). But, unlike in 
the UV regime, b(u) cannot be neglected compared to d(u). In fact, this term is 
crucial for getting a nontrivial solution. In this regime, then, the leading terms in 
equation ( |5.43|) give 

*(«) - 32 ^f^"> (5.50) 

169 (U — Uq) 



which has the solution 

ip (u) « c (w - w ) "^, (5.5i; 



327ri? 6 uJ /2 



where do is an integration constant. Note that the normalization condition remains 
unchanged and cannot be used here because it receives contribution only from the 
UV end due to the exponentially vanishing tachyon potential for large T(u) at the 
IR end. 

Let us now consider the formula, ( |5.38|) , for the lowest mode, using which one 
can compute the eigenvalue Xq. 

/"Umax 1 

W du a(uM(u)) 2 = -Xt (5.52) 

Using a(u) ~ §^1° — , _K in the IR and (|5.51|) in (|5.42|) , we see that rj' (u) ex i/jo(u) 

1/2 

vanishes very rapidly as u — > uq, with a power which grows as u for large u$. 
Moreover, since V(T) vanishes exponentially for large T, the IR region makes a 
negligible contribution to the integral. Therefore, it is reasonable to calculate the 
integral by substituting the UV estimate of the integrand in it. In the UV region, 
o(w) ~ 4 ni5 u A T 2 (u). Moreover, in this region the second term on the right hand 
side of ( |5.42|) can be neglected. So, we get 

1 /"llmax /"Umax fP(ll\ 

-Xt = R e du a(u)( V ' (u)) 2 « #W) 2 / du -U. {^{u)f 

z Juo Ju a\u) 

~ (\^\ 2 f maX 

where Uq > Uq avoids the IR region in the integral and we have defined 



f 2 R 9 

KE W(fl)' (5 ' 53) 

The integral is easily done, giving 

^ _ 1 (T + e- hoUo + T_e feo "°)(T + e- hoUmax + T_e hov ™**) 

f^ g/loC"max — Uq) g — /lo(«max- «o) ' 

From our numerical solutions we see that it is possible to choose uq to be relatively 
large and also satisfy the conditions |T + |e _/,0 "° > \T-\e hoil ° and e ho(llmax ~ fio > > 
e -/i ("max-«o) p or sucn values of the parameters, then, to a good approximation 
( |5.54| ) gives 

A* « -{T + T_ + Tle- 2hou ™*). (5.55) 

Now, let us tune M max to large values. We will do this in a manner consistent 
with the inequality ( 2.19 ). As explained in section O, one way of maintaining this 
inequality is to keep \T + \ and |T_|e /l °" max fixed as u max becomes large. In that case, 



the second term on the right hand side of ( |5.55 ) becomes exponentially smaller than 



the first term as the cut-off is increased beyond some value. We may then neglect 
this term compared with the first term. This gives 



A^ « -T + T_. (5.56) 

Finally, using Aq = m£ and (f4.23|) in this relation, we get 

ml « ^, (5.57) 



This is the well-known Gell-Mann— Oakes— Renner formula, up to a factor of 2. 



6. Summary and Discussion 



This paper further explores our proposal [EI] of a modified SS model, which includes 
the degree of freedom associated with the open string tachyon between the flavour 
branes and antibranes. Here we have extended the analytic treatment of various 
aspects of the problem and supplemented it with extensive numerical calculations. 
We have argued that taking the tachyon into account is essential for the consistency 
of the setup and shown numerically that the solution which includes the tachyon 
is energetically favoured. Our modification preserves the nice geometric picture 
of chiral symmetry breaking of the SS model and at the same time relates chiral 
symmetry breaking to tachyon condensation; the tachyon becomes infinitely large in 
the infrared region where the joining of the flavour branes signals chiral symmetry 
breaking. 

We have identified a parameter in the non-normalizable part of the tachyon 
field profile with the quark mass. It is important to stress that this is the only 
tunable parameter in the modified SS model. It can be traded for the asymptotic 
brane-antibrane separation or the location of the point in the bulk where the brane 
and antibrane join. This provides a natural explanation for the latter parameter, 
which is also present in the SS model, but in that model it doesn't find any 
counterpart in the QCD-like theory at the boundary. In this paper we have 
presented numerical evidence to show that the point where the brane and antibrane 
meet is monotonically shifted towards ultraviolet as we tune the mass parameter 
to larger values. It would seem, therefore, that in our model a brane-antibrane 
pair disappears from the bulk consistently with a quark flavour becoming infinitely 
massive. 

The presence of a non-normalizable part in the tachyon solution requires us to 
introduce an ultraviolet cut-off. The cut-off is needed not only because this part 
grows as one moves towards the ultraviolet region, as in any standard AdS/CFT 
example that includes a non-normalizable solution, but also because the asymptotic 
form of the solution is derived from an approximate equation which is valid only for 



small values of the tachyon. Therefore, the asymptotic solution itself is not valid 
beyond a certain maximum value of the holographic coordinate. We have presented 
sufficient numerical evidence of this phenomenon. Removing the ultraviolet cut-off, 
then, requires tuning the mass parameter to zero. We have explained one scheme by 
which this can be done. This scheme gives an exponential dependence on the cut-off 
to the mass parameter, similar to that discussed recently in |2IJ. The quark mass 
arises from an apparently very different mechanism in this work and the cut-off is 
related to the location of a D6-brane that is present in this model. It would be 
interesting to see if there is any connection between this model and our model. 

Once we have identified the quark mass as a parameter in the non-normalizable 
part of the tachyon, it is natural to expect, by the usual AdS/CFT rules, the 
normalizable part of the tachyon solution to give rise to the chiral condensate. 
To derive an expression for it, however, we need to go beyond the space-time 
independent solutions of section 2. As we have seen, this requires an exact expression 
for the 5-dimensional action for tachyon and brane-antibrane separation fields which 
are now taken to depend on space-time as well as the holographic coordinate. 
We have derived this action in this paper. Using the generalized solutions to the 
equations for this action, then, one can compute the chiral condensate. However, 
one also needs to add counter terms to the boundary brane action to remove from 
it contributions that diverge when the cut-off is removed. 

We have studied in detail the fluctuations of flavour gauge fields on the brane- 
antibrane system. These give rise to vector, axial vector and pseudoscalar towers 
of mesons, which become massive through a kind of higgs mechanism, except for 
the pions. These arise from a gauge-invariant combination of the tachyon phase 
and the longitudinal zero mode of the axial vector field. We have shown that the 
pions remain massless, unless a quark mass (non-normalizable part of the tachyon 
solution) is switched on. For small quark mass, we have derived an expression 
for the mass of the lowest pseudoscalar meson in terms of the chiral condensate 
and shown that it satisfies the Gell-Mann— Oakes— Renner relation. The vector 
and axial vector spectra are expected to be non-degenerate because they arise from 
eigenvalue equations with different tachyon contributions. We have not computed 
these spectra, but it would be interesting to see whether they have the Regge 
behaviour for large masses. 

A non-zero quark mass is essential to correctly reproduce phenomenology in the 
low-energy sector of QCD. Therefore, our modified SS model can be the starting 
point of a more quantitative version of the phenomenology initiated in ]lj. For 
this purpose, our treatment needs to be extended to the non-abelian case, which 
should be a straightforward exercise. The correct tachyon brane-antibrane action 
for curved directions transverse to the branes is not known. It is important to have 
such an action since this would extend the applicability of the present treatment to 
such interesting cases as e.g. the antipodal configuration of the flavour brane system 



and its connection with massless quarks. Another direction in which the present 
ideas can be extended is to discuss this model at finite temperature and describe 
the chiral symmetry restoration transition and study the phase diagram in some 
detail. The connection of chiral symmetry breaking with tachyon condensation 
seems fascinating and a deeper understanding would be useful. Finally, baryons 
have been discussed in the SS model. It turns out that they have a very small size. 
This may change in the presence of the tachyon. This is because in the presence 
of the tachyon, the flavour energy momentum tensor is concentrated away from the 
infrared region where the branes meet. In other words, there is a new scale provided 
by the quark mass. It would be very interesting to investigate whether this effect 
makes any difference to the baryon size. 
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A. Overlapping D8-D8-brane system 

In this case the appropriate DBI action is 

S = - I d 9 a g s V{T) e'* (j-det A L + ^-det A r 

(AU = 9MNd a xfd b x? + K b + \ ((D a r(D b r)* + {D a r) m D b r)) , (A.l) 

where D a r = d a r — 1{Al a — A^ a )r. The classical equation for the profile of the 
magnitude T of the tachyon r can be obtained from ( |2.13| ) by substituting h = in 
it everywhere. We get 

= idJW^, (A . 2) 

where now dx = /(n) _1 w~ 3//2 + T'{u) . In the UV region, assuming T is small for 
large u, we can approximate this equation as 

(VT'(u)y= -7ra 5/2 2», (A.3) 

where we have used the universal small T expansion, V(T) = 7g(l — |T 2 + ••■). 
The general solution 24 to this equation is 

T(u) = u~ 13/s (cicos(4v« /4 ) + c 2 sin(4vW /4 )) + • • • , (A.4) 

24 Equation ( |A.3| ) can be solved exactly in terms of the Bessel functions H^ 1 and H^ . Here we 
give only the leading term. 




where c\ and c 2 are arbitrary constants. Both the independent solutions in this case 
are normalizable, so the interpretation of one of the parameters corresponding to a 
source for the quark mass term is not clear. In view of this, it is not clear how to 



apply the general treatment of |T9| to this case. 

In the IR region, a singular tachyon solution is obtained only for u ~ u k . 
In this region f{u)~ l blows up as (u — Uk)~ l and this drives a singularity in the 
tachyon. Both the potentials in ( [2.5|) and (|2.6|) exhibit singular solutions, although 
the solutions and the nature of singularity are different. For the potential ( |2.5| ) we 
find the solution 

T(u) = U + -—) In- r + 6i + -.-, (A.5) 

V 2 v /m^/ (u - u k ) 

while for (|2.6| ) we get 

T(u) = b 2 (u-u k )- a + ---, (A.6) 

where 61 and b 2 are arbitrary constants and a = 7r ^ Mfc . As in the case with nonzero 
brane-antibrane separation, the IR solution for which the tachyon blows up exhibits 
a smaller number of independent parameters than the UV solution, one in the IR 
as opposed to two in the UV in the present case. A solution with two independent 
parameters in the IR exists (for any potential), but this solution is finite: 

T(u) = T + T 1 (u- u k fl 2 + (-A= + l^(« -«*) + •••■ (A.7) 

3y/u^ 2 V(T ) 

Here T and Ti are the two arbitrary parameters. 

B. Scalar fluctuations 

Here we will assume that T(u,x) and h(u,x) are weakly dependent on x^ and 
expand det(l + K) in Q4.1| ) in powers of space-time derivatives. The action correct 
to quadratic terms in the derivatives is 

S = -f/^/*.»«"V(70V*[l + ^{(l-^)(«T)' 

+ (1 - i£^)l/„3/ W _ £t*L mm }^ (B.l) 



T 



where dx is given by ( |2.11|) , with T{u) replaced by T(u,x) and h{u) by h(u,x). 



Also, the notation (dT) 2 stands for i]' Iu d fl T(u,x)d u T(u,x); similar expressions hold 
for (dh) 2 and (dh).(dT). For the expansion in derivatives to be valid, we must 
require the following conditions to be satisfied: (i) For large values of u, near the 
cut-off -Umax, w e must have \dT\ « ti ma x and \dh\ « 1; (ii) For u ~ u , we must 
have \dT\ « \T'\ ~ (u - u )- 3 and \dh\ « \hT\ ~ (u - n )~ 3/2 . 



Let us now consider small fluctuations around the ^-independent solutions. We 
write T(u,x) = T c (u)+T q (u,x) and h(u,x) = h c (u)+h q (u,x), where T c (u) and h c (u) 
are the x- independent solutions of the classical equations fl2.13| ), ( |2.14| ). We now 
expand the above action and retain only terms up to second order in the fluctuations 
T q (u,x) and h q (u,x). We get 



S 



2V 4 
'1& 



d 4 x 



"0 



V" 



du A\/dr 



VI A 



VI. 



1 + [fT q + j(-fu**h% + T' c T q + h c T 2 h q + h 2 c T c T q )) 



+ {^ + yfi-fu^Kti^ + T' c T' q T q + h c T 2 h q T q + h 2 c T c T 2 ) 



+ ^7^ 3/ X 2 + T ? + T X + hX + ±h c T c h q T q ) - ^(\fu 3/2 KK 



2drA' 



+ T' c T' q + h c T 2 h q + h 2 c T c T q ) 2 + 



u 



-3/2 



2Qo 






rpl2 



in 



3/2 



-h'X(dh q ).(dT q ) + (l- 



\fu^h' 2 l 



-J^ 2 {dh q f 



+ ■■ 



id c c cy q ' v "' v d c 

where we have used the notation V c = V(T C ), d c = dT c , and A = w 13 / 4 V c . As before, 
a prime denotes derivative w.r.t. u, except on V c , for which it denotes a derivative 
w.r.t. its argument. The part of this action linear in fluctuations, Si, which arises 
from the term in the first curly brackets above, is given by 

2V A 



Si 



R 9 



d x I du A 

"0 



^^ c T q+ ^^fu 3 / 2 h' c h' q +T^+h c T 2 h q +h 2 c T c T q ) 
L Vc \/d r 4 



(B.3) 

It is easy to verify that S\ leads to the background equations (|2.13|) and ( |2.14j ). 
This part of the action, therefore, vanishes, except for a boundary term. It is this 
boundary term that gives rise to the chiral condensate. 

The term in the second curly brackets becomes S2, the action quadratic in 
fluctuations, after some manipulations. First, we open the square in the coefficient 
of l/2d 2 term and combine it with the term just before it. That is, we have, 

±-{\ftfl*h* + T' 2 + T 2 h 2 + h\T 2 + Ah c T c h q T q ) 



24 4 



1 
2d r 



- l^{\fu*' 2 ti c ti q + T' c T' q + h c T 2 h q + h 2 c T c T q ) 



1- 



\fu*' 2 h 



/2 



+ 2(2- 



d c 

g£ 2 

dr 



-fu^X + (1 

1 (I 



rp/2 h 2 T 2 

c K + (l-^)(h 2 c T 2 + T 2 h 



'4 J ~ ' q v d ' q ' "' d 



>\) 



'-)h c T c h q T q ) - -{-fu 3 / 2 h' c (Zh' q T q + T 2 h c h q h' q + h 2 c T c h' q T q ) 

+ T' c h c T c (T c T' q h q + h c T' q T q )} (B.4) 

Furthermore, we can rewrite 

<M ) V 2 



A-p=T%T q 

V n V Qj c _ 



vi 



(B.2) 



-A\/dr 



VcK 2 i K (KT C V' T e 

Vc dr. Vc. \ dr. Vr { d t 



r2 



T 2 

-f , (B.5) 



where in the last step we have done an integration by parts over u, used the equation 
of motion ( |2.13 ) for T c , h c and ignored a possible boundary term since it is quadratic 
in fluctuations and so will not contribute to the calculation of the condensate. A 
similar manipulation gives 

A T , h 2 TTT , v f" U "T'\fh 2 c T c \fn\' 



AWd r 



V r . d r . I da il cL 



T 2 



(B.6) 



Combining the above with the other three T 2 /2 terms, we find its net coefficient to 
be 



Aii %-& 



a 



1 



d c 



V'hlT c , hi , V (hlTrS 



'd e + 2-1-2=1 + -*= + 



VcJdr 



'dr 



'dr V d, 



(B.7) 



Similarly, a partial integration using the equation of motion ( |2.14|) allows us to 
combine the two h 2 /2 terms, giving its net coefficient to be 



[l{ h c T^' lfu^h' c | Tl 



d r 



(B.8) 



Collecting all this together, we get the action quadratic in fluctuations: 



S 2 



2V4 



d 4 x / du A 



U() 



l - Cl T 2 + l -c 2 h 2 q + l -c,h' 2 + \c,T' 2 + c 5 h q T q + ceh'X 



1 Q 



C 9 /CVn\2 1 C 10 /OT \ /ryr \ I Cl1 ( XU ^2 



+c 7 h'T q + c 8 h q T' + —^(dT q Y + -^(dh q ).(dT q ) + -^(dh q ) 2 ,(B.9) 



998991 8w 3 Q c ^ 9/ ' 4u 3 Q 
where the coefficients {q} are given by 
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Cl 



rp/2 
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8u 3 Q 
VlhlTr hi TL fhlTrY 



d c + 2-+^== + ^ + ^ c 



°c J -c 



V C J \ d c J ' V c y/dZ \fd~ c ^fd~ c V dc 
h c T 2 \\fu^h' c , T 2 

d r 
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d c 

-d- T A, 

fdr \ 4 ' ' 



/4 v4 
4" 



1 _ 4 J ; ^ )t^ 3/2 , 



KKT* (ri h 2 c T 2 \h c T c 

c 5 = ^^^+ (2 



K vx 



(1. 



/ft 



c 6 



:/: ' V 3/ X, 



4v4 4 



(B.fO) 

(B.H) 
(B.12) 
(B.13) 
(B.14) 
(B.15) 



^4^-^W (B ,6) 



_d c \V c d c J 4 

c 8 = ~J^, (B.17) 

t-i/2 

c 9 " 



4u 3 / 2 V^(l---|-), (B.18) 

c 10 = -u 3 4r h 'X, (B.19) 



/— / 1 -fu 3/2 ti 2 \ 
Pu = ttVv^(l- 4/ \ C ), (B.20) 

with Q c = (1 + fu 3 / 2 h 2 T 2 ). For later convenience, we have explicitly written out a 
factor of l/4w 3 Q c in the coefficients in the last three terms in (P-9| ). 

This action mixes T q and h q and the equations of motion derived from it reflect 
this mixing. After some manipulations, the equations can be cast in the form 

d 2 T q = a{T q + a 2 T' q + a 3 T q + a A h q + a 5 h' q , (B.21) 

d 2 h q = bih q + b 2 ti q + 63^ + b 4 T q + 65T; (B.22) 

where the coefficients {a,j} and {6j}are given by 

a l = Cio(C7 — c 5j + c ll c l) a 2 — Cio(c6 + c 7 — c 8j — C11C4, a 3 = C10C6 — C11C4, 

a 4 = -C10C2 + cn(c 5 - c 8 ), a 5 = C10C3 - cn(c 6 - c 7 + c 8 ), (B.23) 

and 

bl = Cio(Cg — C5) + C9C2, &2 — Cio(c6 — c 7 + c 8j — C 9C3, &3 = C10C6 — C9C3, 

64 = -C10C1 + c 9 (c 5 -c 7 ), 65 = C10C4 -c 9 (c 6 + c 7 -c 8 ). (B.24) 

Here we have used the notation q = (Aci)' /A. As usual, a prime denotes a derivative 
with respect to u. Moreover, d 2 = (—d 2 + <9|) is the flat space-time laplacian. A 
possible term proportional to h!' is not present in ( |B.21| ) because its coefficient, 
(C10C3 — CuCq), vanishes. Similarly, in ( |B.22| ) the term proportional to T" is absent 



because its coefficient, (ci c 4 — c 9 c 6 ), vanishes. 

The equations of motion derived from ( B.9|) are quite complicated in general, 



but they simplify in the two asymptotic regimes of u. 

U — ► U max : In this limit, many of the q are small because they have at least one 
factor of T c or its derivatives in them. The exceptions are c\ ~ /iqM 3//4 , C3 ~ 
ii 9 ' 4 /4, C4 ~ u 3 ' 4 , eg ~ Au 3 ' A and en ~ n 9 ' 4 . Retaining only the dominant terms in 
the equations, we get 

- {u A T q (u, x))' + hi w 4 T q (u, x) = 0, (B.25) 

(u^ti q {u,x)) =0. (B.26) 



The term involving space-time laplacian on the fluctuations can be consistently 
neglected at the leading order since it is non-leading in powers of u, as can be 
verified a posteriori. These equations are identical to ( |2.15| ) and ( |2.17| ) and so have 
solutions similar to ( |2.16| ) and fl2.18| ), but now with parameters that are functions 
of x: 

T q (u,x) = -^(T q+ (x)e- h ° u + T q _(x)e h0U ), 

h q {u, x) = h q0 (x) - h ql (x)u~ 9/2 . (B.27) 

U — > Uq: This limit is more involved, requiring a more detailed analysis. One 
expands T q and h q in powers of e = {u — u Q ) with arbitrary ^-dependent coefficients. 

T q( u , x) = — u^V (^po(x) + e (pi(x) H J , 

h q (u,x) = x p^Uo 3/i e T U (x) + e #i(z) + ••), ( B - 2 8) 

Here, and in the following, we have set /o = 1. One also needs to expand the o^'s 
and frj's in powers of e. Retaining up to the first nonleading power in e, we get 

-3/2 



o. i,, 23e , „.,., 2e N 4m u/ % Vi 23e 

a x =8£e~ 1 (l + -— ), 02 = 2^(1 + — )> «3 = -^~ie\\ + — - 

ILUq Uq TX ILUq 

2ixu n _ 7/9/ _ 65e , AuJ „ i/o,. 21e. 



'; 



a.i 



|H^ 7/2 (1 + ^P), a 5 = ^^ 1 / 2 (l + ^), (B.29) 



and 



26 24m V26 8u 



-i„ 3 ^ , _o^i , 2 ^ , _ W" 2 ^- 236 



6l = _3£ c -i(l+ ), 63 = 2^(1 + — ), b 3 = ^—^(l 



4%" ' ' Uq" TX ' 12Uq 



tx '' 2Auq tx ' 2Auq 

(B.30) 

where £ = — 13mq/8. Substituting these expansions in the equations (|B.21|) , ( B.22| ) 
and comparing different orders of e, we see that a consistent solution exists only for 

us = — 3 and r = —1/2, and then we get 

15 1 Q5Un 

M%) = —7<Po(*)> fi( x ) = £— Vote), &i(x) = — {d 2 + — — )<p {x). (B.31) 
4 6u 8£ 32 

The first of these relations is precisely what is needed to think of the leading terms 
in ( B.28|) as coming from expanding (u — Uq(x)) _1 around a constant uq. The last 
relation shows that when x-dependence is allowed, not all coefficients get uniquely 
determined. In fact, the part of <fo(x) annihilated by the operator on the right hand 
side does not show up in $i(x). 

The above analysis shows that perturbation expansion in "small" fluctuations 
around a constant uq is singular, although we have obtained a solution by a formal 
expansion. 



C. Calculation of the exact (u, x)-dependent action 

This involves calculating the determinant of the matrix (1 + K), whose elements 
are given in (|4.2|) . We will simplify this calculation by making use of the following 
trick. Consider the family of determinants, D(X) = det(l + XK), where A is an 
arbitrary parameter. We actually only need to calculate D(l), but this calculation 
can be reduced essentially to the calculation of the inverse of the matrix (1 + XK), 
which turns out to be much easier than a direct computation of the determinant. 
Consider the following: 

±D(X)=D(XM(1 + XK)- 1 K}. (C.l) 

We can obtain A by integrating this equation, using the boundary condition 
P(0) = 1: 



InD(l) = / dXD(Xy 1 ^-D(X)= f dX tr[(l + XK) 
Jo dX J 



-i 



K] (C.2) 



This reduces the required calculation to finding the inverse matrix M(X) = (1 + 
XK)^ 1 , which may be done as follows. Using the defining equation, (1 + XK)M(X) = 
1, one can express all components of M in terms of M M „: 

M u v = -XK u il M» v , M u u = l-X 2 K\K\M^, M» U = -XM\K\. (C.3) 

Moreover, one can show that M^ v satisfies 

P\M\ = b\, P\ = {8\ + XK", - X 2 K» U K\). (C.4) 

Thus, to find M (A) we need to find the inverse of the P M CT (A) matrix. First note 
that using (|4.2|) we can write 



P'V(A) = *"„ + (3 x {X)d»Td v T + (5 2 {X)d"hd v h + (3 3 {X){d^Td u h + d»hd u T), (C.5) 

where 

a ,^ Xu~V 2 M T'\ Q .^. \f M Ju 3 / 2 h' 2 ^ 

^( A ) = -^7^ 5 P,(X)=p 1 (X)p 2 (X)-(p 3 (X)) 2 . (C.6) 



For A = 1 these reduce to the /3's in ( |4.5|) . Now, from the general structure of the 
P M „ matrix, we can parametrize the M^ u matrix as 

AP„(A) = #*„ + ai {X)d^Td u T + a 2 {X)d^hd u h + a 3 (X)(d»Td u h + d»hd u T). (C.7) 



We have calculated the a's. They work out to be 

1 



cti(A) 
a 2 (A) 
a 3 (A) 



A(A) 

1 
A(A) 

1 
A(A) 



[/3 1 (A)+/3 4 (A)(9/i) 2 ], 
[(3 2 (X) + (3 A (X)(dT) 2 }, 

MX) - p 4 (\)dh.dT]. 



(C.8) 



Here A(A) is a generalization of A defined in ([4.4] ) . It has the same form but with 
the above A-dependent /3's replacing those in (fO]). By definition, A(l) = A. 

Armed with the inverse matrix M(A), we can now compute the trace on the 
right hand side of (|C.2|) . Using (|C.3|) and ( fl2|) , we first note that 



d 



tr[(l + XK^K] = M^(A)— P%(A). 



dA 



(C.9) 



Given the equations (|C.5| )-( |CT8|) , it is straightforward, though tedious, to compute 
the right hand side of the above equation. One gets the simple result 



M^(A)^P%(A) = A(A)- 1 A A (A). 



(CIO) 

UA 

It follows from this and ( |U-2| ) that D(l) = A(l) = A. Hence the complete 5- 
dimensional action is that given in ( |4.3|) . 

To compute the equations of motion for T(u, x) and h(u, x) that follow from this 
action, we will need the following, which can be easily calculated from the relation 
At = dx A and the definition of A given in ( |4.4| ) : 
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(C.ll) 



Using these one can show that 
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We can now give the equations of motion obtained from the action ( |4.3| ): 
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(C.14) 



These can be further simplified using the expressions given in ( |C.11|) , but we will not 
do so here since we will only be interested in a leading solution to these equations in 
the limit u ~ u . As a check, we note that these equations reduce to the equations 
Q2.13D and (2.14) if T and h are x- independent. 
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